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I. INTRODUCTION 

Carrier-phonon interaction is one of the major fac- 
tors that determine the optical and transport proper- 
ties of semiconductors. In semiconductor quantum dots, 
phonon-related effects manifest themselves in a special 
way due to the discrete spectrum of carriers confined in 
these structures. On the one hand, such a discrete den- 
sity of states may slow down carrier relaxation-!. -3 . On 
the other hand, coupling of a confined charge system to 
the nearly dispersionless (thus, also spectrally discrete) 
system of longitudinal optical (LO) phonons leads to the 
formation of correlated carrier-phonon states^ and to a 
reconstruction of the system spectrum, as observed, e.g., 
in intraband absorption experiments^. This effect is a 
generalized form of a polaron, known from bulk systems. 

Apart from the (unmeasurable) energy shifts accompa- 
nying the polaron formation, this effect manifests itself 
by the appearance of polariton-like resonances (anticross- 
ings) between the essentially discrete zero-phonon and 1- 
or 2- phonon states whenever the energy difference be- 
tween two carrier states matches the energy of one or two 
LO phonons^. The 2-phonon resonance is of particular 
interest since the 2-LO-phonon energy of about 70 meV 
falls into the range of typical separations between elec- 
tron energy levels in self-assembled structures. 

Understanding such coherent effects in the carrier- 
phonon coupling in quantum dots is essential not only 
for the correct description of the system spectrum but 
also for the discussion of relaxation properties. In par- 
ticular, the strong coupling between the confined car- 
riers and LO phonons precludes purely LO-phonon- 
mediated relaxation but opens new channels of effi- 
cient two-phonon (acoustic + optical) emission— 9 , which 
may explain efficient carrier relaxation observed in some 
experiments^— . 

Far from the degeneracy (resonance) point, the cou- 
pling to LO phonons may be treated perturbatively. Also 
a first-order polaronic resonance does not present seri- 
ous difficulties, as it involves only a zero-phonon state 
and a one-phonon state. These two states are directly 
coupled by the Frohlich interaction Hamiltonian and the 
resulting resonant anticrossing can be treated, e.g., by a 
second-order Wigner-Brillouin perturbation theor y 13 ! 14 . 



In contrast, the second-order resonance, involving a zero- 
phonon state and a two-phonon state is more complex, 
since the observed anticrossing is due to higher-order, 
indirect couplings via other states. In this case one re- 
sorts to a numerical treatment^&i£~— . The latter is fea- 
sible due to the discrete (dispersionless) character of the 
LO modes, which allows one to describe the LO phonon 
subsystem by a finite set of either orthogonal^ or non- 
orthogonal 1 ^ collective modes. The numerical approach 
to the polaron problem turns out to successfully repro- 
duce the qualitative features. However, the early at- 
tempts to model the experimental results^ suggested 
that certain material parameters must be adjusted in or- 
der to achieve a quantitative agreement with the mea- 
surement data. This seemed plausible 1 ^ as QDs are in- 
herently inhomogeneous systems the quantitative prop- 
erties of which may depend, e.g., on the size and com- 
position characteristics. Such a viewpoint would mean, 
however, that a full understanding and quantitative de- 
scription of the carrier-LO phonon interaction in QDs 
may not be possible. 

The goal of the present paper is twofold. First, we dis- 
cuss the spectrum of the electron-LO-phonon system in 
the vicinity of the two-phonon resonance, including the 
effects of QD asymmetry (ellipticity) , non-parabolicity of 
the confinement potential, as well as the external mag- 
netic field. We show that the experimentally observed 
spectrum of the coupled electron-LO-phonon system can 
be successfully reproduced using only standard material 
constants, without any adjustable parameters. In this 
way we show that the existing theory!^! does allow us 
to understand the physics of carrier-LO phonon interac- 
tions in QDs completely and quantitatively. We calculate 
also the intraband (far infra-red) absorption spectrum of 
a QD ensemble, where the polaronic resonance is clearly 
manifested. Second, we present an effective Hamilto- 
nian approach, based on the quasi-degenerate perturba- 
tion theory which reduces the problem size from a few 
thousands of basis states to just a few and provides some 
insight into the structure of two-phonon polaron states. 
This, in turn, demonstrates that describing the electron- 
LO phonon system does not necessarily have to involve 
heavy numerics and may open the way to the efficient 
description of even higher-order effects. 
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The paper is organized as follows: First, in Sec. UH 
we recall the model of an electron confined in a QD in- 
teracting with LO phonons. Then in Sec. IIII( numerical 
approach is described. The spectrum of a single QD, as 
well as absorption on a QD ensemble and the magne- 
topolaron spectrum are presented in Sec. HV] The two 
following sections fVl and IVT1 describe the effective Hamil- 
tonian approach and present results obtained within this 
method. 



II. THE MODEL 

We consider a self-assembled quantum dot occupied by 
a single electron. The effective confinement potential in 
the xy plane is assumed to be almost axially symmet- 
ric and parabolic, although small corrections accounting 
for ellipticity and non-parabolicity will be taken into ac- 
count. In the z (growth) direction, the confinement is 
assumed to be much stronger, as is typical for these struc- 
tures. The electron confined in the QD is coupled to the 
polarization field associated with the LO phonons. In 
calculations, GaAs parameters will be used. 

The system is described by the Hamiltonian 



The essential part of Hq, accounting for the dynamics 
in the xy plane, is the well known Fock-Darwin Hamil- 
tonian describing a 2-dimensional harmonic oscillator in 
a perpendicular magnetic field. We choose the symmet- 
ric gauge, A = ^(—By,Bx,0), where B is the magnetic 
field, and use the basis of eigenstates of this Hamilto- 
nian (Fock-Darwin states^), denoted as \nm), where 
n = 0, 1, . . . and m = . . . , —1, 0,1,... are the radial and 
angular momentum quantum numbers, respectively (that 
is, hm is the projection of the angular momentum on the 
symmetry axis z). The corresponding wave functions are 



r±\nm) 




' r\m\ 



where (s) is a Laguerre polynomial. Here Ib — 
y^h/(m*uiB) is the in-plane confinement width in the 
magnetic field, where w| = + and oj c — eB/m* 
is the cyclotron frequency in the magnetic field B. In the 
Fock-Darwin basis, the Hamiltonian is 

H = y^^e nm \nm)(nm\, 



H = Hq + H a 



ph ■ 



(1) 



where 



where Hq describes an electron in an isotropic dot, H a 
and H np account for the anisotropy (ellipticity) and non- 
parabolicity of the confinement potential, respectively, 
H ln t describes the electron-phonon coupling and H pn is 
the free LO phonon Hamiltonian. 
The first term in Eq. (TTJ) isIS 



H 



2m' 



(P ~ eA) 



* 2 2 

-m oo r ± 



* 2 2 

-m ui^z 



and describes an electron in an axially symmetric har- 
monic potential in a magnetic field B oriented along 
the symmetry axis, where m* — 0.066m e is the effec- 
tive mass of an electron in GaAs and r± denotes the 
in-plane component of the electron position. The energy 
hujQ is the level spacing for in-plane excitations in the 
absence of magnetic fields. We will refer to this quantity 
as the characteristic energy of the system and use it as 
a system parameter in the discussion that follows. This 
energy parameter is related to the in-plane confinement 
length Zo = y/ h/(m*u>o). The confinement along the z 
axis is much stronger than that in the xy plane, that is, 
uj z 3> ljq. The dynamics along this strongly confined di- 
rection is restricted to the lowest subband, corresponding 
to the ground state wave function 



f 



&V4 



where l z = ^/h/(m*uj z ) is the confinement length in this 
direction. Since the dots we intend to model have the 
height to diameter ratio of about 10 (Ref. 0) we choose 
l z = 0.1/q, which will be fixed throughout the paper. 



e nm = fuv B (2n + \m\ + 1) - -fruj c m. 

The second term in Eq. ([1} describes a weak anisotropy 
(ellipticity) of the confinement potential and has the 
form^S 



ft * 2/ 2 2\ 

-m uj (x - y ) 
2 w B 



^2 v {nm)(n'm')\nm){n'm'\, 



nm.n'm' 



where /3 is a dimensionless parameter and the non- 
vanishing matrix elements in the basis of Fock-Darwin 
states are 



^(02) (00) 



V2 



^(02)(10) — V2, V(oi)(oi) — 1) 



with the symmetries V( 



V, 



nm) (n'm') 



= v, 



(n'm')(nm) 



(nm)(n'm r ) 



Here and throughout the paper, a bar over 
a number denotes a minus sign. This anisotropy term 
leads to the splitting of the p-shell states (in zero mag- 
netic field) given by AE p = (3Hujo which can be read off 
the spectral position of the p-shell states. 

The third part in Eq. ([T]) accounts for non-parabolicity 
of the confining potential 



1 



'0 



where Zo is the in plane confinement length in the absence 
of magnetic field and x *C 1 is a positive parameter defin- 
ing the strength of non-parabolicity. 
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In the basis of Fock-Darwin states, the non- 
parabolicity term reads 

H np = -x^o ^2 V(nm)(n'm')\nm)(ll'm'\, 

(n,rn)(n , m / ) 

where non- vanishing matrix elements are 





V(oo)(oo) 


= 4, 


V(oo)(io) 


= V(io)(o,o) 


= 8, 


V ( oi)(oi) 


= V(oi)(oi) 


= 12, 


V(02)(02) 


= ^(02) (02) 


= 24, 




V(i )(10) 


= 28. 



The electron-phonon coupling Hamiltonian has the 
formal 



-Hint = 



^2 \nm)(n'm'\ ^ F {nm){n , m , } (k)b k + H.c, 



where 



F(nm)(n'm')(k) 

\j 2ve ek I{nm){n m )[q)e 

Here we write the wave vector as A; = 
(fcj_ cos </>, k± sin </>, k z ) and introduce the short-hand 
notation q = /cjJb/2, £, — k z l z /2; f2 is the frequency 
of LO phonons at k — (fiil = 36.7 meV), v is the 
normalization volume for the phonon modes, Eq is the 
vacuum permittivity, and s = (1/Eqo — l/e s ) _1 = 70.3 
is the effective dielectric constant (again, the values 
correspond to GaAs). Note that the Gaussian cut-off at 
k ~ 1/Ib restricts the coupling only to long- wavelength 
modes, so the frequency of LO phonons can be replaced 
by its value at the center of the Brillouin zone (disper- 
sionless approximation). The coupling functions have 
the general symmetry 

F{nm){n'm')(k) = F* n i m i) („ m ) (~ , 

while the / functions for our choice of basis states satisfy 

f(nm)(n'm') n'm')(nm) (q) = f ( («)■ (2) 

The functions f(nm)(n'm')(q) are listed in Tab. |TJ It 
may be interesting to note that the functions F are 
not linearly independent. This follows, e.g., from 
the linear dependence of the subset of / functions 
{/(oo)(oo),/(oi)(oi),/(io)(oo)}, all of which correspond to 
m! — m = 0. The lack of linear independence reduces 
the number of collective modes needed to represent the 
system. 

The last contribution to the Hamiltonian, 

k 

describes free, dispersionless LO phonons. 
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i(q 3 ~ q) 
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10 


-I 2 


i(q 3 - q) 


l-2q 2 + q 4 
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02 


-<?/V2 


iq 3 /V2 


(q 4 - 2g 2 )/ v / 2 


q 4 /2 


02 


-q 2 /V2 


i(q 3 -2q)/V2 (q 4 -2q 2 )/V2 


l-2q 2 + q 4 /2 



TABLE I: Functions /( nm )( n 'm') (?) used in the definition of 
the coupling constants. Functions not listed here can be found 
using Eq. J2). 
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A 


M 
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-q 2 / 


Jxi 


M 


= ±1 


iq/y 


'5a 


M 


= ±2 


-q 2 l 




M 


= ±3 


iq 3 U 




M 


= ±4 


q'/v 





a = B 
— 



a = C 



T 1 



v/5 

i(x 2 q <i -x 4 q) 
y/ x e x 2 —x^x 2 

X4q 4 -X 6 q 2 
\/ x\x % -Xix\ 



(1 - a 2 q 2 + a 4 q 4 )/y/c 



TABLE II: Functions (j>Ma{q) used in the definition of collec- 
tive modes. 



III. THE NUMERICAL APPROACH 

In this Section, we describe the general framework 
for the numerical diagonalization of the carrier-phonon 
Hamiltonian. Then, in Section llVl we present the results 
for a few classes of systems. 

Our approach to the diagonalization of the Hamilto- 
nian ([T]) is based on the collective mode representation of 
the LO phonons^. We use the basis of the electron sub- 
system composed of up to 6 lowest Fock-Darwin states 
(3 lowest energy shells, 2n + |m| + 1 < 3). For this trun- 
cated basis, we define 14 collective phonon modes which 
are needed to exactly represent the carrier-phonon cou- 
pling in the dispersionless approximation, 



Bmc = ^ j;<l>Ma(q)e 



2+lM4> b k , (3) 



where a — A, B,C labels different modes with the same 
angular momentum M and the functions <PMa(q) are 
listed in Tab. [II] For an axially symmetric dot, the ap- 
propriate functions are expressed in terms of the shape- 
dependent parameters (defined for I even) 



1 

47T 3 



^3 q± 

d q^r exp 

q 2 



-> I q\ + \q\ 



4 



Z z /l B =0 



Iz/Ib = 0.1 



16V2TT 



0.0249 
0.0187 



64V2tt 

15 0.0234 



266v^ir 



0.0409 



0.0221 
0.0159 
0.0193 
0.0329 



TABLE III: Comparison between numbers xi calculated for 
strong confinement limit I z /Ib — > and for a realistic value 
Iz/Ib = 0.1 at B = 0. 



Wc define also 



the separation between the unperturbed electron energy 
levels (that is, indirectly, on the QD size). We model 
also the intraband absorption spectrum of an inhomoge- 
neously broadened ensemble of QDs, where the polaron 
resonance is clearly manifested. Next, we calculate the 
spectrum of a single QD with a fixed size as a function of 
an external axial magnetic field and compare the results 
to the existing experimental data. 

When plotting the obtained polaron spectra, we 
present only the states with a sufficient transition proba- 
bility, that way spurious uncoupled states (for the trun- 
cation at the 3-phonon level) are not included. Since we 
focus on the second-order resonance, all important cou- 
plings are contained in the model. 



"2 



a i 



2 



xjx 8 - 
X2X4X6 



2 

X4XQ 
— 



xjx 8 



2 ' 

X/^Xg 



xq — 2a 2 X2 + 204X4 



a\x4 



2a 2 Q4X 6 + a 4 x 6 . 



The numbers xi can be found exactly in the limit of a 
strong vertical confinement, I z /Ib 0, 



X[ 



These limiting values are collected in Tab. IIIII and com- 
pared with those for I z /Ib = 0.1. The leading order cor- 
rection is 0(lz/l B ). 

With the definition ([3]), the collective operators 
Bmu , B\ Ia satisfy the usual bosonic commutation rela- 
tions, [Bmci, B\ Va ,] — Smm'Sccu' (that is, we follow the 
standard approach of orthogonalizcd modest, although 
an alternative approach is also possible 1 -). In terms of 
the collective modes, the interaction Hamiltonian reads 



Mle 2 



21 B £q£ 



nm}(n'm'\ 



X"/(nm)(n'm')a-B m ' — m,a ~\~ H.C., 

where the coupling constants "f( nm )(n'm')a are collected in 
Tab. II VI The mode i?oc couples only to the unit operator 
on the restricted electron subspace, I = J2mn \ fnn)(mn\, 
and can be discarded from the discussion^!. 

The Hamiltonian ([1} is then diagonalized numerically, 
including states with up to 3 phonons, which yields a 
computational basis of 4080 states. The relevance of 4- 
phonon states is discussed in the Appendix. 



IV. NUMERICAL RESULTS 

In this section, we present the results of a numeri- 
cal investigation on the second-order resonant polarons. 
First, we study the system in the absence of a magnetic 
field. We calculate the system spectrum as a function of 



A. Size-dependent polaron spectrum 

We first consider the system without a magnetic field 
and focus on the resonance between the unperturbed ex- 
cited electronic states |0 ± 1) and a set of two-phonon 
states |00,2ph). Here and throughout the paper no- 
tation of this form stands for all two-phonon states 
with the electronic part in the ground state: |00) e i <8> 

(^Ma^M'o' |0) P h). Any other group of states with dif- 
ferent number of phonons is denoted in a similar way. 
The quantities of interest in the present calculation are 
energies of the system eigenstates and the photon ab- 
sorption probability for an intraband transition between 
the ground state and a given excited state. The former 
are obtained directly for the numerical diagonalization, 
while the latter, for an excited state |k), is calculated 
according to 



«i4 _) i*g) 



(4) 



where denotes the numerically calculated ground 
state and d{ is the negative frequency part of the dipole 
moment operator, which depends on the polarization A of 
the optical wave exciting the system. For the a± circular 
polarizations one has (upon truncation to our computa- 
tional space) 



?(-) 



(|01)(00| + V2|02)(01| + |10)(01| 
+V2|01)(02| + |00)(01|)®Iph 



|01)<10| 



and 



(|01)(00| + V2|02)(01| + |10)(01| + |01)(10| 
+\/2|01)(02| + |00)(01|)<g)Iph, 



where I p h is the identity operator on the phonon subsys- 
tem. For the linear polarization along the x and y axes, 
the dipole moment operator is 



oc d ( + } +d[_ \ 



d{: 



otd ( + ] 
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00 A 
B 
C 

01 A 
B 

01 A 
B 
C 

10 A 
B 
C 

02 A 
B 

02 A 
B 
C 



00 

(04X6 — 02X4) / ' \fx~i 
-diyjxg - £Cg/x 4 

-\fx~2 



01 



AE4/2 
AE4/2 



a^XQ / y/XA — («2 + 1)^/3^4 

-anyjxs - x\/xi 

\Jx e - x\/x 2 
\/x s /2 

X 4 /y/2x 2 — \J2x 2 

yJxs/2 - x\/2x 2 



(14 /12 - l)y/X 2 

\]xq- I4/12 

\Jx 6 - I4/12 

(2 — a 2 ) v /ai4 + ( 04 — l)a:6/v / 5?4 
(1 - a 4 )/x 4 yjx 8 - 2jg/a;4 



\/2x 4 — Xq/ y/2x 4 

yj x s /2 - x%l2xi 
\f2x~Z — x&j \J2x 4 
yJxg/2 - x%/2x 4 



s/xi/2 

y/x e /2 
X 4 /y/2x 2 — \J2x 2 

yfxe/2 - x 2 J2x 2 

\J2x 4 — xq/ y/2x 4 
yJx 8 /2 - x\l2xi 

Vx^/2 

(2 — a 2 )yfx4 -f- (cl4 — ^) x q/ \fx~i 
\Jx & - x'l/xii^ - a 4 )/x 4 



TABLE IV: Coupling constants 7( n m)(n'm')a for the collective LO modes. Definitions as in Tab. ITT1 The values not listed in 
the table can be reproduced form the relation 7(„m)( n 'm')a — 7(nm)(n'm')a' 
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FIG. 1: (a,c) Dependence of the unperturbed eigenenergies 
(without electron-phonon coupling) on characteristic energy 
of the QD. (b,d) The QD energy spectrum with electron- 
phonon interaction included. Panels (a) and (b) refer to a 
cylindrically symmetric dot, while (c) and (d) present the 
spectra for an elliptical dot. In (b) and (d) the line thickness 
is proportional to the intraband transition probability |/ifc| 2 , 
the lines are labeled by the polarization of the exciting light. 
Only states for which \fik\ 2 > 0.1 are plotted. 



We investigate two cases: an isotropic QD and an 
anisotropic one (Fig. [T]). In Figs.QJa) andQJc), we have 
presented states from the p shell without phonons as well 
as one-phonon and two-phonon states with the electronic 
part in the ground state (|00, lph), 1 00, 2ph)). The ener- 
gies are shown relative to the system ground state. All 
the other unperturbed states lie in a higher energy range. 
In the case of an isotropic dot, in the absence of the 
electron-phonon coupling, for a certain value of the elec- 



tron level spacing (shaded area in Fig. QJt) , p-shell states 
(|0 ± 1)) intersect with the group of two-phonon states 
with the electron in its ground state (|00, 2ph)). 

When the electron-phonon coupling is included, this 
intersection turns into an avoided crossing pattern, as 
shown for an isotropic dot in Fig. (Ub). Importantly, 
there are no direct matrix elements coupling those states 
and the coupling between the p-shell zero-phonon states 
and the two-phonon states is mediated through one- 
phonon states. Even though this coupling is of the second 
order the anticrossing is quite strong and its width is 1.75 
meV. 

If anisotropy is included, the degeneracy of states 1 01) 
and 1 01) is lifted and two lines in the polaron spectrum 
are observed for different linear polarizations of the inci- 
dent light. As can be seen on Figs. QJc) andQJd), there 
are two intersections [Fig.QJc)] resulting in two anticross- 
ings [Fig. [TJd)] which become visible in the absorption 
spectrum depending on the polarization. 

In the low energy range [36 meV, Figs, [ljb), |T](d)] 
states with small but noticeable transition probability 
can be observed. Their existence is due to the first 
order coupling leading to some transfer of the oscilla- 
tor strength from the p-shell states to 1-phonon states 
[00, lph). 

The position of the resonance in Figs. Hfb) and HJd) 
is shifted with respect to the intersection of decoupled 
states [Figs.[TJa) and[ljc)]. The center of the resonance 
for the interacting case is located at a lower QD energy 
spacing (71 meV) than the intersection point between 
non-interacting states |00,2ph) and |0 ± 1) (73.4 meV). 
Such a behaviour results from the presence of other states 
directly coupled to zero- and two- phonon lines. The 
most important states influencing the position of the res- 
onance are 3-phonon states |00, 3ph) which are relatively 
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close to the 1 00, 2ph) states and effectively reduce their 
energy. 

Since the second-order resonance is relatively strong it 
should also be visible in the absorption spectra of inho- 
mogeneously broadened ensembles. This is discussed in 
Sec. IIVBI More insight to the structure of the second- 
order resonant polarons can be obtained using the effec- 
tive Hamiltonian approach which is presented in Sec. lVIl 



B. Polaron resonance in the ensemble absorption 

In this section, intraband absorption spectra of an in- 
homogencously broadened QD ensemble are calculated. 
QD sizes in self-assembled QD ensembles are always 
given by some distribution. We take this inhomogene- 
ity of sizes into account and theoretically investigate the 
ensemble intraband absorption spectrum in the area of 
the two phonon resonance. 

The QDs are parametrized by their energy spacing 
hojQ l \ where i labels dots in the ensemble. The distri- 
bution of QD energies is assumed to be described by a 
Gaussian function 




40 50 60 70 80 90 100 110 120 
Energy [meV] 

FIG. 2: (Color online) Absorption spectra for two different 
linear polarizations (ax, cry) with anisotropy and without 
non-parabolic corrections. Distribution of characteristic en- 
ergies FkjjQ^ is also presented. The mean energy and its stan- 
dard deviation are taken respectively as ffiUo = 57 meV and 
cr = 6 meV. Labels denote features discussed in the text. The 
dotted vertical line marks the 2-LO phonon energy. The right 
part of the plot is scaled up, as marked. 
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where e = hZJQ and a are the mean transition energy and 
its standard deviation, respectively. Those parameters 
will be chosen in such a way that the second order res- 
onance is located in the high-energy tail of the QD size 
distribution. 

In order to construct high quality absorption spectra, 
we calculate the polaronic states for up to N = 10 5 QDs 

(i) 

with the transition energies Hluq uniformly distributed 
over a sufficiently broad range. From these numerical 
results, the absorption spectrum for the requested polar- 
ization is calculated according to 

N N cv 
i—l k—1 

where Huj is the energy of the absorbed photon, N and 
N cv are the number of QDs used in calculations and 
the number of eigenvalues for i— th QD, respectively, £i K 
stands for K-th eigenenergy of i— th quantum dot and 
|/x K | 2 is given by Eq. (@]). 

The absorption spectrum in the case of a polaron in an 
anisotropic, parabolic confinement potential is presented 
in Fig. [5] The polaronic feature is clearly manifested 
for energies close to the energy of two LO phonons (fea- 
ture B in the plot). Additionally, a phonon replica (C) 
of the main absorption peak is visible for higher ener- 
gies, including the resonance feature (D). The latter is 
slightly broadened and consists of two peaks (D). It is 
worth mentioning that the main absorption feature (A), 
the resonant feature (B), as well as the phonon-replica 
(C) are reproduced correctly in a diagonalization with 
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FIG. 3: (Color online) Absorption spectra on ensembles of 
isotropic QDs for different non-parabolic corrections. Stan- 
dard deviation a = 6 meV, the mean energy fuJo — 62, 67.1 
and 72.8 meV respectively for ascending values of \- 



up to 3-phonon states included. On the other hand, the 
replica of the resonant feature (D) is modelled correctly 
only if 4-phonon states are taken into account (see Ap- 
pendix A). The absorption features for the ax and cry- 
polarizations are not related by symmetry with respect 
to the average QD energy (e). This effect is due to cou- 
pling to a lower lying group of one-phonon states, which 
moves the relevant eigenvalues to higher energies. 

The position of the second-order polaronic feature is 
slightly lower than the energy of two LO phonons 2M7. 
This is mostly due to the interaction with a group of 
3-phonon states |00, 3ph) which effectively reduces the 
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energy of the 2-phonon line. Also the position of the po- 
laronic feature varies slightly with the average value Wjq 
of the QD energy distribution. For the average energies 
higher than the energy of two LO phonons the resonance 
position may shift towards lower energies. However such 
a shift is rather small and does not exceed 1 meV. The 
position of the second-order resonant feature in the en- 
semble absorption is thus not fixed and may vary with 
distribution of QDs in an ensemble. 

Although the potential confining electrons in a QD is 
often considered parabolic, more realistic modelling must 
take into account its non-parabolicity. For the sake of 
simplicity, calculations for a non-parabolic confining po- 
tential are performed for isotropic QDs, that is, the case 
where the anisotropic term H a in Eq. ([T]) is discarded. 
In Fig. 02 we present absorption spectra for different 
strengths of non-parabolicity \- Distribution of the QD 
characteristic energies is tuned so that the main absorp- 
tion peak has the same position. Since in a non-parabolic 
dot higher levels are closer to the resonant group of states 
we expected that the resonance may be shifted down by 
coupling to these states. However, nothing like this is ob- 
served: Even for a strong non-parabolicity, only a negligi- 
ble change in the position of the second-order resonance 
is observed. 



C. Magnetopolaron resonances 

In this section, we investigate a single QD in a mag- 
netic field. We consider a magnetopolaron resonance, 
that is, the case of bringing the state |01) to resonance 
with 2-phonon states |00,2ph) using energy level shifts 
in an external magnetic field. In Figs. EI a ) and Etc), 
we present the spectrum of a single isotropic QD with- 
out electron-phonon coupling in a perpendicular mag- 
netic field, for two different characteristic energies Hu>o. 
In both Figs. Eta) and|4|c) there are several intersections. 
The first one, between the state |01) and a group of 1- 
phonon states, corresponds to the first order resonance. 
States intersecting in the upper part of the charts are the 
purely electronic excited state |01), 2-phonon states with 
the electronic part in the ground state |00, 2ph), and 1- 
phonon states with an excited electronic part 1 01, lph). 
Depending on the size of the QD, those intersections may 
appear at different magnitudes of the magnetic field and 
in different order. 

For the present discussion, the intersection between 
the purely electronic state |01) and a group of 2-phonon 
states 1 00 , 2ph) is of interest [shaded area in Fig.@|a)]. In 
a strong magnetic field, one group of states mediating the 
interaction (|01, lph)) between 0- and 2-phonon states is 
much closer to the resonance than at B = 0. States 
1 01, lph) are directly coupled to both relevant groups 
of states. Since the direct interaction is much stronger 
than indirect one, the resonance between zero-phonon 
and two-phonon states can be strongly intermixed with 
1-phonon states. This is especially the case for QD sizes 
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FIG. 4: (Color online) (a,c) Dependence of the unperturbed 
eigenenergies (without electron-phonon coupling) on mag- 
netic field. (b,d) The QD energy spectrum in magnetic field 
with electron-phonon interaction included. Results obtained 
with numerical diagonalization of full Hamiltonian (lines 
compared with experimental results taken from Refs. H~ 
(points). 



in between those in Figs. Eta) and Etc). 

We compare the results of numerical diagonalization 
with experimental ones 4 ^ obtained for 2 samples differing 
in QDs sizes. The diagrams presented in Figs. Etb) and 
EJd) consist of several lines that represent transitions to 
different excited states. The splitting of the p-shell states 
allows one to uniquely determine the anisotropy param- 
eter (}, while the confinement energy Hloq can be directly 
read off the spectral position of these p-shell states. For 
low magnetic fields, two lines with an opposite Zeeman 
shift emerge. At high magnetic fields, the lower lying line 
strongly interacts with 1-phonon states 1 00 , lph), that is, 
a strong shift towards higher energies is observed. For 
smaller QDs (and thus higher energy spacing Huio), for a 
magnetic field in the range 10-15 T, a clean second-order 
resonance between 1 01 ) and |00, 2ph) states can be ob- 
served [Fig. ED 3 )]. F° r a bigger QD [Fig. H£d)], when the 
upper line comes close to 1- and 2-phonon states, two 
resonances appear. One may assume that the first an- 
ticrossing in the upper part of this chart is due to the 
direct interaction between the states |01,lph) and 1 01) 
whereas the second one is a second-order resonance be- 
tween the states |00, 2ph) and 1 01) . However it must be 
noted that all these states are relatively close to each 
other and both, first and second order polarons can be 
strongly intermixed. 

Previous studie d 4 ' 6 ' 8 ' 18 ' 20 suggested that the electron- 
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Computational 


Resonance 


Resonance 


basis 


width [meV] 


position [T] 


Full model 


o An 
3.00 


IOC 

13. 


3 phonons, 2 shells 


2.10 


12.7 


2 phonons, 3 shells 


2.98 


15.2 


2 phonons, 2 shells 


2.80 


14.9 


Experimental results 


2.9 


12 - 14 



TABLE V: Properties of the resonance area for different trun- 
cations of the computational basis, compared to experimental 
results^. 

phonon interaction, measured by the dimensionless 
Frohlich constant 

e 2 / m* 



needs to be tuned in order to conform with experiment. 
As a result, the Frohlich constant, which in principle de- 
pends only on material properties, was used as an ad- 
justable parameter. Our present results show that such 
a treatment is not necessary. As can be seen in Figs. 
HJb) and0Jd), both resonance positions and widths are 
reproduced without the need to enhance the Frohlich 
constant after a sufficient number of electron shells and 
phonon modes has been included in the computational 
basis. For instance, our numerical solution yields the res- 
onance width of 3.0 meV in the case shown in Fig. QJb), 
which compares quite well with the experimental value 
of 2.9 meV. By performing the diagonalization in re- 
stricted bases, we have found out that the 3-shell, 3- 
phonon model is actually minimal in the sense that a 
further reduction of the basis set leads to incorrect re- 
sults. In particular, leaving out the d-shell states yields a 
strongly underestimated width of the resonance. On the 
other hand, not accounting for 3-phonon states moves the 
resonance to higher magnetic fields and produces nonex- 
istent polaron branches in the resonance area (See Ref.l2ll 
for more details). 

We have also checked to what extent the numerical 
results depend on our choice of the I z /Iq ratio. For I z /Iq — 
0.2, that is, twice larger than used in this paper, only 
small shift in eigenenergies is observed (approximately 
0.1 meV) and the resonance width decreases to 2.8 meV, 
still very close to the experimental result of 2.9 meV. 

V. EFFECTIVE HAMILTONIAN APPROACH 

In the effective Hamiltonian approach, one considers 
the case of an unperturbed Hamiltonian Hq having eigen- 
values E ia grouped into well separated manifolds 2 ^. The 
states can be written in the form \ia), where the Latin 
indices denote different states within a manifold while 
Greek indices refer to different manifolds. The corre- 



sponding energies of the unperturbed system are Ei a . 
Grouping into a manifold means that 

\E ia - E ja \ < \E ia - E jP \ for 

i.e., energy separation between states from different man- 
ifolds, is much larger then energy separations within 
manifold. Moreover, perturbation-induced coupling be- 
tween states from different manifolds should be much 
smaller than the energy separation between these states, 

\(i,a\V\j,fl\ <^\E ia -E j0 \ fora^/3, (5) 

where V is the perturbation. 

The eigenvalue problem is treated by a quasi- 
degenerate perturbation theory. The original Hamilto- 
nian H = Hq + V is transformed into a new one, i? e ff, 
which has no matrix elements between different groups 
of states (up to the required order of approximation) , by 
means of a unitary transformation T = e lS . The ma- 
trix elements of the operator S and of H c g can be found 
iteratively. 

The expression for the second-order effective Hamilto- 
nian has the form 22 

(i\H? s \j) = E ia S ij + (i,a\V\j,a} + 

+ \ E (ha\V\k, 1 )(k, 1 \V\j,a) 



_Ei a Efc^ Eja Ek-y _ 

The first term in Eq. ([6]) represents the unperturbed ener- 
gies, the second one accounts for direct couplings within 
a single manifold, and the last one describes the influence 
of intermediate states from the different manifolds on the 
effective Hamiltonian matrix. This last term represents 
indirect second-order couplings between the states of the 
manifold of interest which result from the couplings to 
other manifolds eliminated by the unitary transforma- 
tion T. 

The effective Hamiltonian method is a powerful tool 
for calculations and interpretation of various systems. In 
the present case of a two-phonon polaron resonance, it 
is very helpful since in the resonance area a purely elec- 
tronic state and a group of two-phonon states form a well 
separated manifold. The method can be used for a de- 
scription of both magnetopolaron resonances and a size- 
dependent polaron spectrum, though it is more accurate 
in the latter case, since the condition for the appropriate 
relations between energy spacings and coupling strengths 
[Eq. <j5j) ] is fulfilled in this case with a greater precision. 

If one considers the case of the size dependent spec- 
trum (without a magnetic field) the energy difference 
between states from different manifolds at the point of 
the resonance is at least Kl = 36.7 meV. On the other 
hand, if the purely electronic state 1 01) is brought to res- 
onance with the 2-phonon line |00, 2ph) using a magnetic 
field, the energy separation is significantly lower. For the 
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cases presented in Figs. |4}a), EJc), it is approximately 15 
and 30% smaller, respectively (assuming that the coupled 
state 1 01, lph) is considered as a member of the mani- 
fold). If we take into consideration that the average di- 
rect coupling between states from different manifolds is 
about 3 meV, even 30% reduction in energy separation 
might impair the applicability condition. 

The effective Hamiltonian approach automatically in- 
cludes the nested coupling polaron structure^, since 
from the whole spectrum of the states it filters out the 
ones that are coupled to each other through intermediate 
states (for a given truncation of the basis). 

VI. RESULTS: EFFECTIVE HAMILTONIAN 

In the following section, we investigate the second- 
order resonance using the effective Hamiltonian ap- 
proach. We consider the case of size dependent spectra 
as well as a magnetopolaron resonance. 

We apply an appropriate treatment to obtain the ef- 
fective Hamiltonian in both cases, although the latter 
one might be less applicable in the case of larger QDs 
(Zeeman tuning brings one-phonon states close to two- 
phonon line intermixing first and second order resonances 
and thus breaks the condition for sufficient manifolds sep- 
aration). Since the energy separation between different 
states depends on the size of the QD each case should be 
studied individually. 

The effective Hamiltonian approach allows us to gain 
some information about the second-order polaron struc- 
ture. Since we choose relevant indirectly coupled states 
only we can get much insight into the mediated interac- 
tion between 0- and 2- phonon states. As the effective 
Hamiltonian matrix is much smaller than the full Hamil- 
tonian matrix its diagonalization is much faster and in- 
terpretation of the spectra is easier. Contributions from 
different intermediate states can easily be separated and 
studied. In particular, we investigate the influence of d- 
shell and 3-phonon states on the polaron spectra and on 
the coupling strengths appearing in the effective Hamil- 
tonian. We show that the quasi-degenerate perturbation 
theory not only allows one to describe the resonance area 
in detail, but also explains why both d-shell and 3-phonon 
states have to be used in order to correctly model second- 
order polarons. 



A. Polaron resonance at B = 

In this section, we consider the effective Hamiltonian 
approach to the second-order resonance between states 
1 01} and 1 00, 2ph) in the absence of a magnetic field. For 
the sake of simplicity, we assume that the QD is isotropic 
(it is always possible to introduce anisotropy perturba- 
tively) . 

Although, for our truncated basis, there are 105 two- 
phonon states only 6 of them couple indirectly (in the sec- 



matrix elements 


intermediate states 


s 


P 


d 


lph 


3ph 


(OlltfefflOO,^ 1 ) 





82.8% 


17.2% 


1 





(Oll^efflOO,^ 1 ) 


88 


.9% 


11.1% 


1 





(OlItfefflOO^B 1 } 








100% 


1 





(OlltfefflOO,^ 3 ) 








100% 


1 





<01|tf cff |00, A 2 B l ) 








100% 


1 





(01\H cH \00,B°A 1 ) 


50% 


-50% 





1 





(00,A°A 1 \H e g\00,B°A 1 ) 


100% 








0.5 


-0.5 



TABLE VI: Relative influence of different electronic shells and 
n-phonon states on the effective Hamiltonian couplings in the 
limit Iz/Ib — > 0. In the second line, the contributions from 
the two shells are taken together as they partly cancel each 
other. 



H e B matrix elements 


coupling strength 


(OlItfefflOO^A 1 ) 

(OlltfefflOO,^ 1 ) 
(OllffefflOO,^ 1 ) 

(01|F e ff|00,A 3 A 5 ) 
(01 H cff 00, A 2 B l ) 


29v / 3/(768v / 2T) ~ 0.0261 
13v / 3/(384v / 2^ : ) « 0.0234 

-1/(128Vtt) « -0.0044 
-5 v / 5/(256 x /2li : ) « -0.0035 

-1/(256VtF) « -0.0022 



TABLE VII: Indirect coupling strengths in the effective 
Hamiltonian approach. Coupling strength expressed in units 
of e 2 /{2le e) = 32.7 meV. 



ond order approximation) to the purely electronic state 
1 01) . There is a small number (12) of intermediate states 
which produce nonzero couplings in the effective Hamil- 
tonian. Their contribution is presented (grouped by shell 
and number of phonons) in Tab. lVIl If the numerical val- 
ues for different shells or different phonon numbers sum 
up to zero the relevant states are decoupled. 

As we can see in Tab. IVI1 taking the d-shell into 
account not only introduces three additional couplings 
(rows 3 to 5), but also changes the strength of existing 
ones (increase of over 17% and 11% in rows 1 and 2, re- 
spectively). What is more important, taking 3-phonon 
states into account completely decouples the 2-phonon 
state |00, B a A 1 ) as shown in row 7. Here, A 1 denotes 
the phonon mode created by the collective operator B[ A 
[Eq. ©], etc. Although Tab. ED is constructed for the 
limiting case of an infinitely flat QD (l z — > 0) and at ex- 
act resonance (ftwo — 2M7) this decoupling is preserved 
also for non-flat QDs and is an effect of equal spacing 
of the electronic eigenenergies in a parabolic confining 
potential. 

The properties of indirect couplings via 1-phonon and 
3-phonon states, revealed by the effective Hamiltonian 
structure, allow one to qualitatively understand why a 
3-shell, 3-phonon model is required for correct modelling 
of the resonant polaron spectrum. Leaving out 3-phonon 
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FIG. 5: (Color online) Comparison between effective Hamil- 
tonian approach and full Hamiltonian diagonalization. Influ- 
ence of disregarding 3-phonon states on resonance area also 
included. 



matrix elements 


coupling strength 


d— shell influence 


(OlltfefflOO,^ 1 ) 


1.786 


7% 


(OlliJefflOO,^ 1 ) 


0.627 


11% 


(0l\H cS \00,A 3 A*) 


-0.157 


100% 


(OlItfrfflOO^B 1 ) 


-0.143 


100% 


(01\H cB \00,A 2 B 1 ) 


-0.056 


100% 



TABLE VIII: Indirect effective Hamiltonian couplings for case 
of magnetic field tuning. Values obtained for electron exci- 
tation energy huo = 60 meV in a magnetic field of 14.1 T, 
expressed in meV. 



approach are nearly the same as those obtained with full 
Hamiltonian diagonalization. The largest shift between 
eigenenergies found with those two methods in the range 
of the QD energy spacing presented in Fig. [5] is lower 
than 0.16 meV. 



states leads to a reconstruction of the spectrum (addi- 
tional lines in the theoretically modelled absorption spec- 
trum are observed, see Fig. [5]) since the coupling be- 
tween states \00,B°A 1 ) and |00, A A 1 ) in the absence of 
3-phonon states is relatively strong. For the case of a flat 
QD, the coupling strength factor is 0.0288 which is the 
strongest coupling comparing to other ones shown in ta- 
ble lVIII However, if the 3-phonon states are included this 
coupling vanishes completely. This shows that a trunca- 
tion of the computational basis can not only cancel some 
indirect couplings^ but can also lead to the appearance 
of nonexistent ones. This increases the number of opti- 
cally active states in the resonance area and affects the 
spectrum both qualitatively and quantitatively. 

Although the cancellation of an indirect interaction 
between certain 2-phonon states in the presence of 3- 
phonon states is quite general and similar cancellation of 
the mediated interaction between other n-phonon states 
(n > 0) can be observed it does not necessarily trans- 
late to the reduction of the size of the relevant polaronic 
subspace. In general, such a decoupling is a quantitative 
effect and may depend not only on the structure of the 
model, but also on the values of the couplings. For that 
reason the quasi-degenerate perturbation theory seems 
to be the right approach, which takes in to account the 
nested coupling structure^, as well as system dependent 
quantitative effects. 

The influence of rf-shcll and 3-phonons states on the 
diagonal elements of the effective Hamiltonian is found 
to be less than 0.4% and 2.5% respectively. Coupling 
with the two-phonon state |00, A 2 B 1 ) is very weak and 
it is possible to discard it from consideration. Such a 
reduction of the computational basis does not produce 
any important effects. The area of the second-order res- 
onant polarons can now be described by a Hamiltonian 
of dimension 5. 

The results obtained with the effective Hamiltonian 



B. Magnetopolaron resonance 

The effective Hamiltonian treatment for a non-zero 
magnetic field differs slightly from the previous case of 
zero magnetic field, since the paramagnetic term brings 
the one-phonon state |0l,A 2 ) close to the second-order 
resonance. This additional state couples to the purely 
electronic state and to certain two-phonon states and, 
depending on the QD size, can be of major importance. 
Nevertheless, the resonance area can be precisely de- 
scribed using the quasi-degenerated perturbation theory 
even in the case of relatively big QDs. When the pre- 
viously mentioned one-phonon state is too close to the 
energy of two LO phonons it simply needs to be included 
as a member of the considered manifold. Since both sam- 
ples [Figs. H{b), HJd)] consist of relatively large QDs this 
one-phonon state has to be included in the group of rel- 
evant states. 

In the construction of the effective Hamiltonian, a 
set of 6 indirectly coupled states {|01), lOO,^ ^ 1 }, 
lOO,^ 1 ), \00, A 3 A 2 ), lOO,^ 1 ), lOO,^ 2 ^ 1 )} is chosen 
as a basis (decoupling of certain 2-phonon states, dis- 
cussed in the previous section, is already taken into ac- 
count). At this point, we are interested in the strengths 
of indirect coupling mediated by 1-phonon states, so we 
temporarily exclude the state |01, A 2 ) from the manifold. 
In this way, the influence of all intermediate 1-phonon 
states on the second-order resonance can easily be deter- 
mined. To keep the model reasonably accurate in the ab- 
sence of this state Tab. I Villi is calculated for a QD with 
a slightly larger characteristic energy Hu>o, so that the 
energy separation between the relevant manifold and the 
state |01,A 2 ) is also larger. Indirect coupling strengths 
for the resonance condition with the influence of the d- 
shell are presented in Tab. I Villi As one can see, in the 
case of the magnetopolaron spectrum, the presence of 
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FIG. 6: (Color online) Polaron eigenenergies in the resonance 
area. Dots: experimental results from Refs. lines: results 
obtained with the effective Hamiltonian approach. Samples 
and QD material parameters as in Fig. 2] 



FIG. 7: (Color online) Absorption spectra on ensemble of 
isotropic QDs obtained with and without 4-phonon states. 
Standard deviation 6 meV, the mean value KuJo = 62 meV. 



the d-shell increases the coupling strength significantly 
(rows 1-2) and couples 3 additional 2-phonon states with 
the state |01) (rows 3-5). The character of the influence 
of 3-phonon states is the same as that discussed in the 
previous subsection. 

As we pointed out earlier, in the case of relatively 
large QDs (when the characteristic energy is lower than 
60 meV), additional one-phonon state in the effective 
Hamiltonian basis needs to be included. If this is done 
then both methods: diagonalization of the full Hamil- 
tonian and quasi-degenerated perturbation theory ap- 
proach produce nearly the same results in the resonance 
area with differences in the obtained eigenenergies and 
resonance widths smaller than 0.15 meV and 0.01 meV, 
respectively. Comparison between the results obtained 
with the effective Hamiltonian approach and the exper- 
imental ones are presented in Fig. [5] Good agreement 
between the theory and experiment can be observed. 



VII. CONCLUSIONS 

In this paper, we have theoretically studied the reso- 
nant features in the spectrum of an electron confined in 
a self- assembled QD and interacting with LO phonons. 
We have focused on the second-order resonance induced 
by the indirect interaction between the first excited elec- 
tronic shell (p-shell) and the electronic ground state with 
two LO phonons. We have studied this second order 
resonant polaron spectrum as a function of the dot size 
(energy level separation) and external magnetic field. We 
have also calculated the absorption spectra for an inho- 
mogeneous ensemble of QDs and shown that polaronic 
feature is clearly manifested in these spectra. Our re- 
sults, compared to the existing experimental data, show 
that a properly constructed model is able to quantita- 
tively reproduce the observed polaron resonance without 
any need for free or adjustable parameters describing the 
interaction between confined electrons and LO phonons, 
except for shape and size parameters that can uniquely 



be extracted from the intraband absorption spectrum. 

In order to get more insight into the structure of the 
polaron spectrum in the resonance area, we have devel- 
oped an effective Hamiltonian approach based on a quasi- 
degenerate perturbation theory. We have shown that by 
tracing the structure of indirect couplings mediated by 
1- and 3-phonon states, a very small set of relevant basis 
states can be identified which span the space of resonant 
polaron states. 

The presented results show that the spectrum of the 
coupled electron-LO phonon system can be reliably mod- 
elled based on the standard theories and computational 
techniques developed for confined systems. Moreover, 
they demonstrate that this modeling may be consider- 
ably simplified by applying perturbation theory methods, 
without loosing the accuracy of the results. 

Appendix A: Influence of 4-phonon states on the 
ensemble absorption 

Since the 4-phonon states do not couple 0- and 2- 
phonon states, their influence on the 2-phonon feature 
is negligible. On the other hand, they are important 
if one considers one-phonon replica of the second-order 
polarons. We compare here the absorption spectra for 
two computational bases: one including only states with 
up to 3 phonons and the other one with additional 4- 
phonon states (Fig. [7]). For the sake of simplicity, results 
were obtained for the case of an isotropic QD without 
anharmonicity corrections. 

The influence of 4-phonon states on the main absorp- 
tion peak (A), as well as on 2-phonon resonance (B), is 
marginal. On the other hand, 4-phonon states are of 
major importance for one phonon replica features (C,D). 
Since those features appear as a result of the interac- 
tion between 1-phonon and 3-phonon states, 4-phonon 
states have similar influence on them as 3-phonon states 
had on the interaction between 0- and 2- phonon states 
(decoupling of previously strongly coupled states). As 
a result of introducing additional interacting states, a 
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quantitative change is observed in the intensity of the 
phonon replica (C) of the main absorption peak. The 
most important change in the ensemble absorption is re- 
lated to the phonon replica of the resonant feature (D). 
If 4-phonon states are omitted this feature is broader and 
consists of two peaks. On the other hand, if we include 
4-phonon states, the replica of the resonant feature con- 



sists of one sharp peak, which is shifted towards lower 
energies. The energy shift is mostly due to presence of 
directly coupled 4-phonon states with energy AMI which 
is higher than the energy of the feature D, whereas the 
change in its shape is related to a reconstruction of the 
spectrum in the area of the feature D, introduced by those 
additional states. 
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